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Abstract 

A closed system of equations for the local Bloch vectors and spin correlation functions of three 
magnetic qudits, which are in an arbitrary, time-dependent, external magnetic field, is obtained 
using decomplexification of the Liouville-von Neumann equation. The algorithm of the derivation of 
the dynamic equations is presented. In the basis convenient for the important physical applications 
structure constants of algebra su(2S-|-l) are calculated. 
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I. INTRODUCTION 



At present it is not known which experiment will lead to the reliable, prototypical quan- 
tum computing device. Quantum systems with two states, called qubits, are taken to be the 
basic unit for quantum information processing. The quantum systems realised on coupled 
qubits are widely investigated for the purpose of creation of a quantum computer. In the 
same time there is a possibility to solve this problem, using the systems consisting from cou- 
pled qudits (multilevel systems). Qudits possess a number of properties which differ from 
properties qubits and may have advantages for quantum information processing. For exam- 
ple, two qutrits can be entangled more strongly than two qubits The main purpose of 
this study is the derivation of dynamic equations for three qudit system, which is placed in a 
magnetic field, for entanglement study [2] and other applications. The article is organized as 
follows. Section II contains the model Hamiltonian including the three-particle interactions. 
In section III the Liouville-von Neumann equation for the matrix density of three qudits in 
the variable magnetic field is obtained in the Bloch's representation in terms of local vectors 
and spin correlation functions taking into account three-particle interactions. We describe 
conservation laws which control numerical calculations effectively. The Appendix describes 
the Hermitian basis and analytic formulas for structure constants. 



II. MODEL HAMILTONIAN 



The Hamiltonian of three coupled different qudits (particles with spin 1/2, 1, 3/2, 2, ...) 
in the external ac magnetic field h = {hi, h2, h^), where /ii,/i2,/i3 are the Cartesian com- 
ponents of the external magnetic field, we will present in a form of decomposition on a 
complete set orthogonal Hermitian matrices Ca, operating in the Hilbert space of three 
qudits H^^ ® H^^ ® H^^ with {2Si + 1)(252 + 1)(2S3 + 1) dimension 

^ = ^ VrCf ® ® Cf^ (1) 

where Si, S2, S3, are qudit spins, ® - denotes a direct product. C^^, Cf^ are the Hermi- 
tian matrices operating in spaces H''^^, H^^ , H^-^ accordingly (see Appendix). The functions 
hgcTT contain one- , two- and three-qudit interactions 
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III. DECOMPLEXIFICATIN OF THE LIOUVILLE-VON NEUMANN EQUA- 
TION 



The Liouville-von Neumann equation for the density matrix p, describing the dynamics 
of a 3-qudit system, has the form 

idtp = [H,pi p(0) = po, p+ = p, Trp = 1. (2) 

It is convenient to rewrite the equation Eq. ([2]) having presented the density matrix p 
as well as the Hamiltonian H, in the form of decomposition on a complete set orthogonal 
Hermitian matrices C„ 



where Ci^2,3 = y '^^'^'^'''^g'^'^"'"'^'' . We define the three coherence Bloch vectors -Rmoo, -Rono, -Roop 
, which are widely used in the theory of magnetic resonance characterize the local properties 
of the individual qudits, 

ciC2C3i?„oo = Tr pC^,^ ® C^' ® (4a) 

ciC2CsRono = Tr p C^^ ® ^ ® Cq^^ , (4b) 
CiC2CsRoop = Tr p Cq^^ ® C^' ® . (4c) 

while the tensors RmnO, RmOp, Ronp, Rmnp 

ciC2C3Rmno = Tr p C^' ® ® Co^^ , (5a) 

CiC^CsRmOp = Tr p C|,^ ® C^' ® , (5b) 
ciC2C3Ronp = Tr p Co^^ ® ^ ® , (5c) 
CiC2C3i?^„p = Tr p C^^ ® ® , (5d) 



describe the spin correlations. In the formulas Eq. (jlj), Eq. (jS]) Cq' = Cq ^ = y ^^"-^ — ^-^Si) 
is the unit matrix in dimension {2Si + 1) and each of the Latin indices designates a 
set of matrices C^l s- ,C^i. c. ,Cfj. , where 1 < A;*^' < 25^, and 1 < q^^ < k^' in 
C'fc^,,^,^,, Cls,^^s,,y and 1 < fc^' < 2^„ in Cf,, in steps of 1, ^ = 1, 2, 3. 

Let's formulate the linear algorithm of a decomplexification of the Liouville-von Neumann 



equation. 

1. Insert the Hamiltonian Eq. ([1]) and the density matrix Eq. ([3]), decomposed on Hermitian 
basis into the Liouville-von Neumann equation Eq. ([2]). 

2. Multiply the Liouville-von Neumann equation by all elements of the basis C^^ ®Cg^ ®C^^ 
in turn. 

3. Execute operation of a trace taking for each equation. 

4. Apply the formula Tr (Cf ^ ® ® C^'){C^' ® Cf^ ® C-f ® ® ) = 

r^5i (^Si (^52^52/^-52 /^Ss/^Sa^Sa 

II II Ly p II Ly^ . 

5. Express a trace from the three matrices through structure constants according to formu- 
las Eqs. (IA4yA5l) . 

6. As a result the use of the structure constants symmetry the real terms in each equation 
are cancelled out, and the purely imaginary terms are duplicated. 

7. The imaginary unit i is cancelled. 

This algorithm is easy to apply for the system of more than 3 qudits. 

In terms of the functions Rap-y the Liouville-von Neumann equation becomes real and com- 
prises a closed system of first-order differential equations for the local Bloch vectors and 
spin correlation functions 

dtRmoo = C2C3ejim{hjooRm + hjkoRiko + hjokRiok + hjkrRikr), (6a) 
dtRono = CiC3ef^^{hoiQRQko + hji^Rj^o + hQijRQkj + hjirRjkr), (6b) 
dtRoop = CiC2efjp{hf)(jiRQQj + hkoiRkoj + hQkiRokj + hikiRikj), (6c) 



dfRmnO = C2C3ejf^hjQQRinQ + CiC^C -IJlQjQRmkO + C2C3^jim^jnoRm + ClC3Ggkn^'rnqoRokO + 
C3(e^m5'fcgn + 9jim(^kln)hjqoRikO + 02036 jl^hjOgRinq + OiC3e^lJlojqRmkq + 
+ 02036 jl^hjnqRiOq + OiC3e^^,JlmiqRokq + C3 (Cj/^^'j^^^ + 9 jimOlkn)^jk^ikq, 

(6d) 

dtRmOp = 02036 jl^hjQoRiQp + OiC26^^phoQkRmOq + O2O36jl^hj0pRi0Q + 01626^^^ JlrnOkRoOq'^ 

02{6jlrn9kqp + 9jimOkqp)hj0kRi0q + 02036jl^hjkoRikp + OiC26fg^phoklRmkq+ (6e) 
Ol026f^phmklRokq + 02036jl^hjkpRikO + C2(e^im5'fgp + 9jim^rqp)^jkrRikq, 
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dtRonp — ClCs^ikn^OioRokp + ClC2e^p/loOi-Roni3 + ClCsC^^^/loip-RofcO + CiC2e^^phonkRoOq + 

ci(ef4^Jp + gfkn^j'JhoijRokg + ciCsel'^haoRikp + ciC2eg,/iioii?mg+ (6f) 

S2 ( S'2 S S2 S \ 

cic^CglnhigpRiko + ciC2eig'phiniRi0q + Ci{ei^^g,.^p + gij^^e.^^p)hiriRikq, 

dtRmnp = C2C3^jimhjnpRm + ClCs^gkn^rnqpRokO + ClC2^kqp^'rnnkRoOq + C2C3^^im^j0pRin0 + 
ClC^e^jlJlQjpRmkO + C?.{<i^jlm9lkn + 9^jim^%n)^jQpRikQ + C2C^e^jl^hjnoRiOp + 
ClC2ellphonkRmOq + C2 (c'^'^Sf;^^ + gflm'^fgp)hjnlRiOq + CiC3e^^^hmioRokp+ 
CiC2ef^phm0iR0nq + Cl(ef45'zgp + gikn^^qp)hmilRokq + C2Cze^jlmh jOoRinp + 
ClCse^ln^OjoRmkp + Czi^^jlmQlkn + djim'^fkJ^jloRikp + CiC2efgJlQQjRmnq + 
C2{ejljn9lqp + 9jim'^lqp)hj0lRinq + Cl{^%n9lqp + 9jkn'^lqp)hojlRmkq + 
i^jim9lkn9rqp ~^ 9jim^lkn9rqp ~^ 9jim9lkn^rqp ~ ^jim^lkn^rqp)^jlrRikq- 

(6g) 

As idtp^ = [H, p"] {n = 1,2,3,...) at unitary evolution there is the numerable number of 



motion, 



fl. 



From 



conservation laws Trp = Ci = 1, Trp^ = C2, . . . , where C„ are the constants o 
from which only the first (25*1 + 1)(25'2 + 1)(2S'3 + 1) are algebraically independent 

dcf 

the conservation of purity, for which (p^)ifc = ip)ik, the polynomial (square-law) invariants 
are obtained. The square polynomials also control the signs Rap-y- In the external dc field 
the energy of system is the constant: 

E = Tr Hp. (7) 
Unitary evolution preserves the length of the generalized Bloch vector Jj^^^^^Ss 



JJS1S2S3 _ . Jl'^^^ -)- + RqQp + -RmnO + R^nOp + R'onp + Rmnp- (8) 



The qudit with {2Si + 1) states in the environment of two other qudits is described by the 
reduced matrix p^^ 

P^^ ^ — Toe — TTT^-^oooCo'^ + -RmooCm (9) 
in which i?ooo = 1, and functions -Rmoo are determinated by the system solution Eq. ([6]), as 
the equations for the reduced matrices are not closed. 

For two different coupled qudits we have H = ^ha/sC^^ ® C^^, 



The dynamic equation Eq. ([2]) for two different qudits takes on the real form in terms of the 
functions Rmo, Rom, Rmn SiS a closed system of differential equations Eqs. flll|12lll3l) for the 
set of initial conditions: 

dtRmo = C2epl^{hpoRio + hpiRu), (11) 
dtRom = Ciepf^{hopRoi + hipRu), (12) 

dtRmn = ^pim [^2 (/'-pn-RjO + hp^Ri^) + g^j^^hprRnj + 

^pin [ci{hmpRoi + hopRmi) + grlm^rpRli\ , (13) 

where by definition 



TipC^^^Cl' =C,C2Rat3. (14) 

The functions Rmo, Rom describe the individual qudits and the functions Rmn define their 
correlations. The Liouville-von Neumann equation for one qudit takes on the real form in 



terms of the functions Rj as a closed system of differential equations 



dtRi = e%hiRj. (15) 



The set of equations for 3 qubits has been obtained in [8l 

The set of equations Eq. (I6l) with the initial conditions has wide applications, since 
the magnetic field enters in the form of arbitrary functions. It allows to make numerical 
calculations for continuous (a paramagnetic resonance in a continuous mode), as well as 
for pulse modes (a nuclear magnetic resonance). By means of this system it is possible to 
investigate the entanglement dynamics of qudits in a magnetic field as the entanglement 
measures are expressed in terms of the reduced density matrices or of populations. Another 



important application of the system Eq. ([6]) is quantum approach to the Carnot cycle 
10|. 11], [l2|, when a working body is a finite spin chain. 



IV. CONCLUSION 



The simple algorithm of the derivation of equation system for coupled qudits, which 
are in an arbitrary, time-dependent external magnetic field has been presented. It is not 
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necessary for the basis to be Hermitian since the results of calculations are independent of 
the choice of base, but there is the main advantage with the Hermitian basis. It is that 
the Liouville-von Neuman equation not involve any complex numbers and can be solved 
using real algebra. This is not true for non-Hermitian bases. Real algebra makes numerical 
calculations faster and simplifies the interpretation of the equation system Eqs . (jH]). 



This basis forms a natural basis for calculations on coupled spin systems 14j because all 
the single-spin operators are part of the complete basis when the unit operator is part of 
the single-spin basis. 
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APPENDIX A 

Let {Cf , C|', } be a base of su(2S+l) algebra, where S = 1/2, 1, 3/2, ... is the spin 
quantum number, n = {2S + 1)^ — 1. We have according to [l^ 

Cf C/ = + zf^.Cl Tr Cf = 0, Tr C^Cf = c5,,, (Al) 



-'j d ijk^k ' -^^ ^' ^ j 

^ijk ~ 9ijk + '^^ij 



^ijk — 9ijk + '^^ijk-: (^2) 



hence 



z[Cf , C/] = 2eJ,Cf , {Cf , C/} = -^E5., + 24.,Cf , (A3) 
1 

2ic 



4^ = ^Tr[Cf,C/]Cf, (A4) 



4-. = ^Tr{Cf,Cf}Cf, (A5) 

where d = 2S +1, E the unit matrix in dimension x c?, c is a constant, Tr is a symbol for 
trace. It is easy to see that the structure constants e^-^, and g^^^. are completely antisymmetric 
and symmetric in the displacement of any pair of indices. 
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1. Hermitian basis 



The structure constants of su{2S + 1) algebra have important physical applications. In 
order to calculate the structure constants we have to choose the basis. The basis is based on 
linear combinations of irreducible tensor operators. The matrix representations of irreducible 
tensor operators T^^ [3] can be calculated using the Wigner 3jm symbols: 

s 

T,^^= v/(25 + l)(2fc + l) i-lf~"'{'J,i')\S,m><S,ml (A6) 

m,m'=—S 

where < A; < 25*, and —k < q < k in steps of 1. The normalization is such that T^q = E. 
It is known that the Cartesian product operators Sx, Sy, and 5*2 for spin S = ^ are Hermitian 



and can be calculated from irreducible tensor operators 17l | 

1 

2^2 



Sl = -^(Titi -Th = ^r \, (A7a) 



1 






2 


(: 











H 


i 





Sy' = ^(Tl^i + ^I'l) = ^ r / h (A7b) 



1 11 1 /I , 

Si = 7T^i!o = 7;\ • (A7c) 



2 2 



\0 -1 



Allard and Hard 13] have formed linear combinations of the irreducible tensor operators 
not only for single- quantum coherences, but also for all coherences according to 



Cl,x = \/ ^^^^^(^t, + i-inij, q + 0, (A8a) 



Cl,y = nh^^^iT,',-, - i-^YT,',), Q ^ 0, (A8b) 



Cl. = \ ' , ' T,',o, g = 0, > 1, (A8c) 



Ct = \l ' ; ' E, (A9) 



_ . IS{S+ 1] 

where 1 < A; < 25, and 1 < g < A; in Cig^,Clgy and 1 < A; < 2S, in C^^ in steps of 1. 
The matrices Eqs. (IAS I) are traceless and their number is equal to (25* +1)^ — 1. Using the 
well-known relations for the irreducible tensor operators 

(TiX = i-^yTk,-,, (AlO) 



we can see that matrices Eqs. (lASP are Hermitian. Using the formula from 16 1 



TTT'T,',. = {-iy(2S+l)S,,k'5, 



q,-q 



(All) 



it is easy to show that the basis is normahzed so that = C'f^., Sy = Cf , Sz = Cf 



irrespective of the spin quantum number S, i.e. 



1)(2S+ 1) 



(A12) 



The set Eqs. (lA8fA9p is complete. The matrices C^^z are diagonal 



(A13) 



There also exist the other useful bases [l8 , 19|. From the physical point of view, for 
important physical applications the basis [14!] is preferred. 



2. Analytic formulas for structure constants 



There are 27 combinations in threes including the repetitions: XX'X", XX'Y" 
XX' Z"..., where X = C^^^, X' = Cf,^,,., Y" = Z" = Cf,,, and so on. 

use of the symmetrical properties of the Wigner 3jm symbols and the formula 2.4 (23) 



rT\. rp S rpS rpS _( -| \ 25+fe + fc' + fc'7 9 Q , 1 

-■-J- ^ k,q-'- k',q'''' k",q" — \ ^) "r ^ 



[(2A: + l)(2fc' + l)(2A:"+l)]Mr;r} 



k k' k"^ I k k' k" 
q q' q" 



(AM) 



where {| |.' |"} is the 6j symbol, allows us after substitution of Eqs. flASP in Eq. (j4]), Eq. ([5]) 
, to calculate all structure constants of su{2S + 1) algebra. Let us introduce the function 

F{k, k', k", S) = ^^^S{S + 1)(25 + 1)(2A; + l){2k' + 1)(2A;" + 1) {| f }. (A15) 
v3 

All antisymmetric structure constants are zero for K = k + k' + k" even and nonvanishing 
antisymmetric structure constants in terms of 3jm and 6j symbols have the explicit form 
are presented by formulas Eqs. (]A16a|A16b|A16cl) for K odd: 



,5 

^XX'Y" 



F 

71 



(-1)^ i^q-V-::) + i-iy' {'q + (-1)^" 



k k' k" 
q q'-q" 



(Al6a) 



72 



k k' k' 
-q q' q" 



k k' k' 
q -q' q" 



k k' k' 
q q'-q'' 



(Al6b) 
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e'xy'z" = -F{-iy(iy,). (Al6c) 

All symmetric structure constants are zero for K = k + k' + k" odd and nonvanishing 
symmetric structure constants in terms of 3jm and 6j symbols have the explicit form are 
presented by formulas Eqs. (]A17allA17bllA17c|) for K even: 



9xx'x' 



9XY'Y" 



F r 

7! 



V2 



( T \ q ( k k' k' 



k k' k' 



(Al7a) 
(Al7b) 



q") + (-1) 

X'Z" = a'yy^z" = n-^y (; 4' o") > ^^Iz'Z" = ^(-1)" (o o'o") • (Al7c) 

We have in X, F 1 < A;, k' , k" < 2S, 1 < q < k,l < q' < k', 1 < q" < k" and in Z 
1 < k,k',k" < 2S in steps of 1. 

The straightforward calculation confirms that the structure constants efjj^ and gfji^ are com- 
pletely antisymmetric and symmetric in the displacement of any pair of operators. In other 
words it is Cxx'y" = ~'=^xy"X'^ 9xx'Z" = 9yy'Z" so on. 
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